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Abstract: These lecture notes focus on the bound state sector of QCD. Motivated by 
data which suggests that the strong coupling a s (Q) freezes at low Q, and by similarities 
between the spectra of hadrons and atoms, I discuss if and how QCD bound states may 
be treated perturbatively. I recall the basic principles of perturbative gauge theory bound 
states at lowest order in the h expansion. Born level amplitudes are insensitive to the ie 
prescription of propagators, which allows to eliminate the Z-diagrams of relativistic, time- 
ordered Coulomb interactions. The Dirac wave function thus describes a single electron 
which propagates forward in time only, even though the bound state has any number 
of pair constituents when Feynman propagators are used. In the absence of an external 
potential, states that are bound by the Coulomb attraction of their constituents can be 
analogously described using only their valence degrees of freedom. The instantaneous 
A field is determined by Gauss' law for each wave function component, i.e., for each 
position of the valence constituents. Solutions for A obtained with a boundary condition 
that imposes an asymptotically constant energy density give rise to a linear potential for 
color singlet qq and qqq states. The strength of the linear potential is determined by the 
boundary condition and is of lower order in a s than the gluon exchange interaction, which 
may then be treated as a higher order perturbative correction. Bound states evaluated 
to a given order in a s and h must have the full symmetry of the exact theory, including 
the dynamic boost invariance of states quantized at equal time. The wave functions are 
indeed found to have such a hidden invariance, which ensures the correct dependence of 
the energy eigenvalues on the center of mass momentum. Thus relativistic bound states 
can be studied using perturbative methods. 
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1. The strong interactions 

There is a broad consensus that Quantum Chromodynamics is the theory underlying the 
strong interactions. This conviction first arose from general features of QCD which agree 
with observations. The realization that the QCD coupling a s decreases with momentum 
transfer (asymptotic freedom) was decisive for establishing QCD as a serious candidate 
theory in 1973. QCD was also seen to have approximate isospin and chiral symmetry for 
small quark masses. The theory later has had an impressive success in describing a large 
variety of data on short-distance processes. These involve small a s and thus directly probe 
the QCD lagrangian via the perturbative expansion. Numerical methods based on discrete 
lattice approximations of space-time have demonstrated that the properties of QCD at 
long distances, in particular its bound states (the hadrons), are at least qualitatively in 
agreement with observations. 

What I summarize above in a few sentences would provide material for many lecture 
courses. The vast amount of information has fortunately inspired a "Resource Letter" 
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[1], which at 39 pages is no letter but a comprehensive collection of references to the QCD 
literature. It covers all aspects of QCD and guides to material organized by scope and depth 
via brief remarks. Many excellent introductory lectures on QCD have been published as 
well, see for example [2, 3, 4]. I shall not try to rival these here. Instead, I have chosen to 
focus on a less well covered but fascinating aspect: the bound states. 

Color confinement and the spontaneous breaking of chiral symmetry are novel features 
of QCD, whose origin is still poorly understood. Both are closely connected to the physics 
of hadrons as bound states of QCD. Color confinement implies that only color neutral 
states (hadrons and nuclei) propagate macroscopic distances and thus reach experimental 
detectors. The spontaneous breaking of chiral symmetry explains the absence of parity 
doubling in the hadron spectrum and gives a reason for the small mass of the pion. Hadrons 
also have other puzzling properties: Hard scattering shows that they have a large (infinite) 
number of highly relativistic quark and gluon constituents, whereas the hadron spectrum 
reflects just the valence quark (qq or qqq) degrees of freedom [5]. This property makes 
hadrons quite different from more familiar bound states, such as atoms and molecules, 
whose spectra reflect the relative motion (vibrations, rotations,...) of all their constituents. 

The u, d valence quark masses contribute < 2% to the nucleon mass. Hadrons are in 
fact the only truly relativistic bound states found in Nature. In QED atoms the relativistic 
effects are of O (a) and can be calculated perturbatively with exquisite precision [6]. A 
relativistic description of hadrons would seem to require a s > 1 and thus methods beyond 
perturbation theory. Together with the increase of the running coupling a s (Q) at low 
momentum scales Q ("infrared slavery") this has led to the general view that QCD is 
strongly coupled at distances ~ 1 fm. 

But data suggests otherwise. I review in Sec. 2 the indications that a s (Q) freezes at 
a modest value ao ~ 0.5 as Q — > and that QCD remains perturbative at long distance. 
There is little freedom in a perturbative approach since the lagrangian determines the 
expansion nearly uniquely, and qualitatively correct results should be obtained already 
at lowest order in a s . Hence the relevance of perturbation theory for soft QCD physics 
is an assertion that can be addressed - and potentially excluded. Nevertheless, aspects 
of relativistic bound states that definitely can be addressed perturbatively receive little 
attention: 

- Does the Hydrogen atom wave function Lorentz contract? [7] 

- How does the Dirac wave function describe higher (e + e~ pair) Fock states? [8] 

- Is there a Born term for bound states, as there is for scattering amplitudes? [9] 

Modern courses in field theory devote surprisingly little attention to bound states. The 
perturbative description of bound states differs qualitatively from that of scattering ampli- 
tudes, since bound state poles arise through the divergence of the perturbative expansion. 
It thus seems worthwhile to devote these lectures to bound state issues. I shall not at- 
tempt an objective overview, but rather focus on issues which seem relevant for a possible 
description of QCD bound states. 
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Figure 1: Left: The dependence of the QCD coupling a s on the momentum scale Q [10]. Right: 
The average coupling ao defined by (2.1) as determined from data in the range < Q < 2 GeV 
[13]. 

I discuss the evidence that QCD remains perturbative at long distances in Sec. 2. 
Sec. 3 introduces general issues related to the perturbative description of bound states. The 
possibility of using retarded propagators at Born level and how this affects wave functions 
is discussed in Sec. 4, and then applied to Dirac bound states in Sec. 5. Sec. 6 addresses 
the dynamical (as opposed to manifest) boost invariance of states defined at equal time of 
their constituents. Each order in an h expansion must be Lorentz invariant, as discussed 
in Sec. 7. In gauge theories the ^4° potential is at each instant of time determined by 
the charges through Gauss' law. In Sec. 8 I note that the Coulomb field of the Hydrogen 
atom therefore is different for each wave function component (position of the charges). 
This allows me in Sec. 9 to consider solving Gauss' law with a boundary condition where 
A does not vanish at spatial infinity. This generates a linear confining potential without 
breaking the symmetries of the theory. In Sees. 10 and 11 I use this boundary condition to 
derive the QCD bound state equations at lowest order in h and a s for qq and qqq states. 
Final comments are given in Sec. 12. 

2. The freezing of a s (Q) 

Strong interaction phenomena are characterized by the scale Aqqd — 200 MeV ~ 1 fm _1 . 
This fundamental constant also determines the value of a s in the perturbative regime: 
a s (M z ) = 0.1184 ± 0.0007 [10]. The running of a s (Q) shown in Fig. 1 (left panel) agrees 
well with data down to the scale of the r lepton mass, where a s (1.8 GeV) ~ 0.33. 

It is often claimed that a s (Q) grows large as Q — > Aqcd, and that confinement is 
a consequence of strongly coupled QCD. Actually, several analyses indicate that a s (Q) 
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"freezes", i.e., becomes independent of Q at low scales [11]. A dispersive approach [12] 
which uses moments of event shapes to extract an average coupling ceo at low scales, 

1 

«o(w) = — / dQa s (Q) (2.1) 
fJ-i Jo 

gave [13] ao — 0.5 (for fij = 2 GeV) as shown in the right panel of Fig. 1. A recent 
analysis of event shapes in e + e~ annihilations which combined the dispersive method with 
perturbative calculations at NNLO accuracy resulted in [14] 

a s (M z ) = 0.1153 ± 0.0017 (exp) ± 0.0023 (th) 

a = 0.5132 ± 0.0115 (exp) ± 0.0381 (th) (2.2) 

A moderate coupling at long distances would help explain why perturbative results are 
found to qualitatively describe data at low Q 2 . Examples include Bloom-Gilman duality 
[15], precocious dimensional scaling [16] and the inclusive distribution of hadrons in and 
between jets [17]. Hadrons are, like atoms, classified according to the spin and orbital 
angular momentum of their valence constituents. The success of the quark model in de- 
scribing hadron masses, magnetic moments and other properties suggests a perturbative 
structure. In the words of a well-known expert [18]: 

QCD is about to undergo a faith transition. 

QCD practitioners prepare themselves - slowly but steadily - to start using, 
in earnest, the language of quarks and gluons down into the region of small 
characteristic momenta - "large distances". 

This is a perplexing statement since we know that QCD is not perturbative in the same 
sense as QED is - there are no free quarks or gluons. Applied to the phenomenological 
quark model potential 

4 rv 

V{r) = kr - — (2.3) 
3 r 

it suggests that the coefficient k of the linear term is of lower order in a s than the gluon 
exchange contribution — |a s /r. In Sec. 9 I show that a linear term arises in the pertur- 
bative expansion if a non-vanishing energy density oc Aq CD is imposed on the solution of 
Gauss' law. Then k oc y/cts Aq CD is indeed of lower order than the O (a s ) gluon exchange 
contribution, which may be treated as a higher order perturbative correction. 

3. Bound states in perturbation theory 

Bound states appear as poles in scattering amplitudes. In a perturbative expansion the 
poles are generated by the divergence of the sum - no finite order Feynman diagram has 
a bound state singularity. This is intuitively understandable, since the constituents must 
continuously interact in order to stay bound. It may nevertheless seem surprising that 
the QED perturbative series diverges even for atoms, however small is a (I return to 
this in Sec. 7). Having to sum an infinite set of diagrams raises the question of which 
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diagrams to include in the sum: Unlike in ordinary perturbation theory one cannot simply 
order the calculation according to the power of the coupling a. Summing different sets 
of diagrams leads to different approximations for the bound states, which a priori appear 
equally justified [19]. 

The perturbative expansion of scattering amplitudes can be viewed not only as a 
power series in a, but also as a loop expansion. Loops are associated with powers of the 
Planck constant h [20]. Born terms may have different powers of a depending on the 
process, but are always of lowest order in h in the sense that they have no loops. The very 
successful QED calculations of non-relativistic atoms are also based on a loop expansion. 
The loops are counted in the interaction kernel, which is then iterated as a geometric series 
to generate the bound states. Does this mean that the Schrodinger equation, which follows 
from iterating single photon exchange with no loop correction, is the Born approximation 
for non-relativistic bound states? And if so, can this concept of a Born term be extended 
to relativistic bound states, defining a specific lowest order approximation (bound state 
equation)? The h expansion has received rather little attention in a field theory context, 
but it appears that the answer is positive to both questions. I discuss the h expansion 
further in Sec. 7. 

Having identified the Born term one may ask under what circumstances this approx- 
imation to a bound state is reliable. Unsurprisingly, it turns out that this requires the 
coupling a to be small enough for a perturbative (loop) expansion of the kernel to be well 
motivated. At first this appears to make the Born term concept useless for relativistic 
states. In QED atoms the constituents move with speeds v/c ~ a. Hence relativistic 
motion requires a large coupling a > 1, which invalidates the loop expansion. The Dirac 
equation with a —a/r Coulomb potential in fact gives complex bound state energies for 
a > 1 [20]. 

A systematic study of bound states using perturbation theory requires expanding both 
in h (which at lowest order defines the Born term) and in the coupling a (or a s in QCD). 
The general understanding (which I review in Sec. 7) is that these two expansions are 
coupled: Each loop adds a power of h and a power of a. At small a we then have 
non-relativistic states described by the Schrodinger equation - unless there is a confining 
potential at lowest order in a. In Sec. 9 I argue that the perturbative expansion does allow 
a linear potential, through a non-vanishing boundary condition in Gauss' law. 

In gauge theories the ^4° field does not propagate (there is no 8qA° term in the la- 
grangian, unless it is introduced via the gauge condition). Therefore A is determined by 
the positions of the charges, separately for each component of the wave function and at 
each instant of time. Imposing a non- vanishing boundary condition at spatial infinity when 
solving for ^4° gives rise to the linear term in the potential (2.3), with k determined by the 
asymptotic field strength. Thus the QCD scale Aqcd may appear in the solution of the 
equations of motion even though the equations themselves have no such parameter. 

The possibility of studying the properties of relativistic bound states within the very 
constrained framework of perturbation theory is interesting in its own right. Bound states 
of lowest order in a and h retain the Lorentz and gauge invariance of the theory. The 
dynamical boost invariance of the equal-time bound states found using this method is non- 
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Figure 2: The tree diagram in 4-momentum (p°,p) space on the left splits into two time-ordered 
diagrams in (t,p) space (right hand side) when using the Feynman propagator (4.2). With the 
retarded propagator (4.3) both the positive and negative energy contributions are given by the first 
diagram on the right, and the Z-diagram is absent. 

trivial and (to my knowledge) unique. In considering the physical relevance of such an 
approach it is worth recalling that we rarely reflect on the reason (in fact we do not know) 
why QED does not confine electric charge. This ignorance has not kept us from doing 
perturbative calculations, implicitly assuming A to vanish at large distances. The validity 
of this expansion for QED is certified by its precise agreement with data. 

4. Retarded vs. Feynman propagation 

Tree (or Born) diagrams usually give a good first approximation to scattering amplitudes. 
They are special in the sense that all internal propagators are off-shell and hence insensitive 
to the ie prescription of the propagators. Specifically, Feynman and retarded fermion 
propagators 

n f + m e 

Sf/r{p ' p) = ! FvWTW (41) 

give the same result for Born terms, whereas Feynman propagators are required in loop 
integrals. The choice of propagator nevertheless makes a difference even at Born level when 
the interactions are time-ordered (p° — > t). The E < components propagate backward in 
time when using the Feynman prescription, 

S F (t, V) = i^ [0(t)(E pl ° - p ■ 7 + m e )e^ - 0(-t)(E pl ° + p-j- m e )e lE ^\ (4.2) 

which gives rise to Z-diagrams as shown in Fig. 2. With the retarded prescription there is 
only forward propagation, 

S R (t, P ) = ^ [(£ P 7° - V ■ 7 + m e )e- iE ^ + (E p7 ° +p- 1 - m e )e lE ^] (4.3) 

The first diagram on the rhs. of Fig. 2 thus gives the entire time-ordered contribution when 
the retarded prescription is used, and the Z-contribution vanishes. For either prescription 
the sum of all time-ordered diagrams (integrated over the interaction times) gives the same 
Born amplitude. 

There is usually no need to time-order the interactions of scattering amplitudes, hence 
also the issue of ie prescription is irrelevant at Born level. Time-ordering is, however, 
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Figure 3: Electron scattering from a static external potential. The p° component of the electron 
momentum does not change during the scattering. The initial and final electron momenta are 
denoted (p°,0) and (p°,p), respectively. 

required to describe bound states using wave functions that give the spatial distribution of 
the constituents at an instant of time. When using Feynman propagators the infinite set of 
diagrams that generates a bound state pole has contributions with an unlimited number of 
nested Z-diagrams. This means that the wave function has components with any number 
of constituents. 

As an illustration consider an electron in a static Coulomb field (Fig. 3). The energy 
component p° of the electron's 4-momentum is unchanged by the interactions, which only 
transfer 3-momentum. In the absence of loop corrections the electron's Green function is 
given by the sum 



where S is the free propagator and K the Coulomb interaction kernel — ie^/ A°(k) (the 
product involves a convolution over k). The last equality in (4.4) displays a bound state 
pole contribution. It is easy to verify that the residue R(Er,p) satisfies the Dirac equation 
in momentum space, i.e., we have Dirac bound states. 

For p° > the electron propagators S in each diagram of Fig. 3 are insensitive to 
the ie prescription at the negative energy pole since p° ^ —E p in (4.1). Consequently the 
positions p° = Er of the bound state poles in G(p°,p) are exactly the same whether we use 
Sf or Sr propagators. Yet if we time-order the interaction vertices these two propagators 
will generate very different diagrams. Each Feynman propagator gives a Z-contribution 
as in Fig. 2, while with the retarded propagator only the single electron (of positive or 
negative energy) is present at any intermediate time. The equal-time wave functions of the 
bound state are thus quite different in the two cases. The wave function we usually solve 
using the Dirac equation has only the single electron degree of freedom and corresponds 
to the use of a retarded electron propagator. The wave function of Dirac bound states 
corresponding to the Feynman propagator has an unlimited number of Fock components 
and is seldom (ever?) given explicitly. 

5. The Dirac wave function 

The bound states \tp) of non-relativistic quantum mechanics are usually determined as 
eigenstates of the hamiltonian, 



G(p°, p) = S + SKS + SKSKS + . . . = S + SKG 



R{E R ,p) 
P°~E r 



+ ... 



(4.4) 



H\tp)=E\<p) 




The same condition holds in relativistic field theory, but is impractical even for the simple 
case of an electron bound by a Coulomb potential just discussed. Since the hamiltonian can 
create an e + e~ pair from the Coulomb field the state \(p) must necessarily contain an infinite 
number of pairs - the same conclusion that we reached previously. Hence in relativistic 
theory one usually determines bound states as poles of Green functions in 4-momentum 
space rather than using a hamiltonian formulation. 

According to our discussion above the standard (single particle) Dirac wave function 
follows from using retarded propagators in a perturbative evaluation of Green functions 
as in (4.4). We may then ask whether we can set the boundary conditions in the oper- 
ator equation (5.1) correspondingly, such that it defines a single particle bound state \tp) 
with a Dirac wave function. For this we need a "retarded vacuum" in which the electron 
propagator 

S R (x -y) = fl (0| T[i/>(x)${y)] \0) R (5.2) 
agrees with the retarded one in (4.3). The definition 

|0)« = iV- i n<AlO) (5-3) 

where all the positron states are filled 1 and iV is an (infinite) normalization constant works 
since 

Hx)\0) R = (5.4) 

implies no contribution for x° < y° in (5.2). A single-electron state with both positive and 
negative energy components can then be parametrized by a Dirac (c-number spinor) wave 
function ip{x) as 

\<p,t) = J d 3 xiji{t,x)v a (x)\0) R (5.5) 

where a sum over the Dirac index a is implied. With the QED hamiltonian in the Inter- 
action Picture 

H(t) = j d 3 x ${t, x) [ - iV • 7 + m + ej°A°(x)] ip(t, x) (5.6) 

the state (5.5) in (5.1) gives, using {if) a (t, x), tpg(t, a;')} = 5 3 (x — x')8 a p and (5.4), the 
Dirac equation for the wave function <p(x) of a bound state of energy E in the external 
Coulomb potential A°(x), 

(-iV • 7 + ej°A°(x) + m)(p(x) = Ej°(p(x) (5.7) 

In Sec. 4 we saw that the energy E is independent of the boundary condition (Feynman 
or retarded) when the potential is static and there are no loop corrections. The possibility 
to describe the same state \(p) using different wave functions is not so surprising when we 
recall that the time-ordering of events which are separated by a space-like distance depends 

1 Equivalently. in the retarded vacuum all the negative energy states are empty, which is why pair 
production is suppressed. 
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on the observer. Thus the relative magnitude of the two time-ordered diagrams in Fig. 2 
depends on the frame, only their sum is Lorentz-invariant 2 . In the infinite momentum 
frame the Z-diagram vanishes altogether and one arrives at a limiting picture where the 
hamiltonian does not create particles from the perturbative vacuum. Here I propose to 
eliminate the Z-diagrams in a different way, by using retarded fermion propagators at Born 
level. The bound states are found to be Lorentz-covariant and can thus be boosted to an 
arbitrary frame. To this lowest order contribution perturbative corrections may then be 
systematically applied. 

6. Dynamical boost invariance 

Many approaches to relativistic bound states strive to retain manifest Lorentz covariance. 
A two-particle Fock state wave function for a bound state \(f,P) with CM momentum P 
is covariantly defined by 

$5j(x - y ) = (o\T{Mx)My)}\v, P) (6.1) 

This Bethe-Salpeter wave function has a 'kinematical' frame dependence, 

q> p '( x > - y >) = S(A)$ p (x - y)S~ 1 (A) (6.2) 

where S(A) is the usual Dirac matrix representation of the Lorentz transformation A, 
P' = AP and x' — y' = A(x — y). Eq. (6.2) relates a Bethe-Salpeter wave function defined 
at equal time (x° = y°) to wave functions which generally are at unequal time ix' ^ y'°)- 
Equal-time wave functions are related to each other in a non-trivial way since the notion 
of equal time is frame-dependent. 

Quantum field theory is generally not manifestly covariant since one quantizes on a 
space-like surface, usually the one of equal time [21]. Poincare invariance is ensured by 
requiring the transformation generators to satisfy the corresponding Lie algebra. Gener- 
ators which preserve the quantization plane (such as space translations and rotations for 
equal-time surfaces) define 'kinematic' transformations for which the symmetry is man- 
ifest. Generators which transform the quantization plane (such as time translation and 
boosts for equal-time surfaces) involve interactions and correspond to 'dynamical', hidden 
symmetries. 

A distinguishing feature of bound states is that they have a stationary time depen- 
dence, i.e., they are eigenstates of the hamiltonian as in (5.1). Their wave functions define 
the spatial distribution of the constituents at an instant of time. Boosting equal-time states 
is a dynamical transformation which is generally as difficult as diagonalizing the hamilto- 
nian. This and the problem of unlimited particle production discussed above challenge 
hamiltonian approaches to relativistic bound states. 

Classical relativity would suggest that equal-time wave functions Lorentz contract like 
a rigid rod, which is indeed often claimed without further discussion. This issue was only 

2 Provided that also the external potential is Lorentz-transformed. Lorentz covariance is more properly a 
feature of states that are bound by the mutual interactions of their constituents and experience no external 
force. I consider this case in the next Section. 
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addressed recently in the basic case of the Hydrogen atom [7]. It turned out that the 
wave function of the moving atom receives a contribution from both the e~p and the e~p7 
Fock states. In the atomic rest frame the instantaneous Coulomb field ^4° dominates and 
only the electron and proton constituents are present. When these charged constituents 
move relativistically with the CM they generate a propagating vector field A implying 
intermediate states containing (longitudinal) photons. The e~p7 component of the atomic 
wave function is required for the energy eigenvalue of the Hydrogen atom to have the 
correct dependence on its CM momentum, E = \J P 2 + M 2 , but does not Lorentz contract 
according to classical relativity. 

In the previous section we saw how a relativistic electron bound by a static Coulomb 
potential is described by the single particle Dirac wave function when one uses retarded 
propagators. The bound state energies are insensitive to the choice of propagator when 
only tree diagrams are retained, which we referred to as the Born approximation. Provided 
this approximation indeed defines the lowest order in an h expansion it must share all the 
symmetries, and in particular dynamical boost invariance, with the exact result. Hence we 
turn now to a discussion of the h expansion in field theory. 



7. The h expansion 



It is generally accepted that each loop in a Feynman diagram brings a factor of h [20, 22], 
and that physics in the limit of h — > is classical. However, the first statement has been 
challenged [23] , and the second statement also needs to be stated more precisely. Classical 
quantities such as mass and charge appear, for dimensional reasons, multiplied by powers 
of h in the lagrangian. In order to properly define an k — > limit one must therefore specify 
the behaviour of all parameters in the theory [9]. 

It is instructive to consider the H — > limit of the harmonic oscillator in quantum 
mechanics. The propagation of a particle from to (tf,Xf) is given by the path 

integral 



A(xi,xf,tf -ti) = / [T>x(t)]exp 



[£>£(*)] exp 



im 




2h 




im 




IT. 


It, 



dt(x 2 - urV) 



dt{e 



(7.1) 



In the second expression all explicit dependence on h was removed by a mere redefinition of 
the variables, £ = xj \fh. Hence it is clear that the full quantum structure of the harmonic 
oscillator, including its bound states, remains in the h — > limit. The usual argument for 
an approach to classical physics is that the rapid variation of the phase factor exp(iS/h) 
selects the classical path for which the action S is stationary. This constraint is evaded for 
paths whose length scales as x oc since then S oc h. The bound state wave functions are 
in this category since their extent scales with £. Classical physics dominates in the h — > 
limit only for paths whose length is independent of h. 

In field theory one usually sets h = 1 at the outset, causing the dimensions of length 
L and energy E to be inversely related. For the present discussion [9] we wish to keep the 
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Figure 4: The first two ladder diagrams contributing to non-relativistic atoms in the limit of small 
coupling a. 



factors of h explicit, with dimension [h] = E ■ L. We keep c = 1, so the units of space and 
time are both L. The QED action 



~>QED 



(7.2) 



should have the same dimension as h, hence [ip] = E l l 2 L~ l and [A] = i? 1 / 2 L~ 1 / 2 . From 
the fact that the fine structure constant a = e 2 /(4-7r/i) is dimensionless we conclude that 
the dimension of the classical charge is [e] = E^-^L 1 / 2 . This is why the charge in the QED 
action (7.2) enters in the form e/h. Similarly the classical mass appears as m/H. 

In considering the classical limit of quantum field theory it seems natural to keep the 
classical quantities e and m fixed as h — > [23]. For a connection between the power of h 
and the number of loops we need instead to keep e = e/h and rh = m/h fixed [20]. After 
a rescaling of the fields, ip = i^/Vh and A = A/y/h we get 



i/j — eVhjji — mj i/j (7.3) 



where h appears exclusively in the combination ey/h. Adding a loop to any given Feynman 
diagram gives a factor (e\/h) 2 and thus a factor of h. This establishes the equivalence 
between the h and loop expansions 3 . The conclusion is the same for QCD. 

How about the h expansion for bound states? In the weak coupling limit (a — > 0) 
QED atoms are given by a sum of ladder diagrams, the first two of which are shown in 
Fig. 4. Each ladder increases the number of loops by one, so atoms would seem to be of 
no definite order in h. On the other hand, the binding energies E n = —^a 2 m e /n 2 are of 
definite power of a = e 2 h/4-7r. The loop and h expansions appear not to be equivalent in 
this case. 

Bound states are indeed a special case because the momenta in the ladders depend 
on a. For the ladder sum to maintain an overlap with the bound state wave function the 
momentum transfer in Fig. 4 must be of the order of the Bohr momentum, \q\ = O (am e ), 
while the energy transfer should be commensurate with the binding energy, q° = O (ct 2 m e ) . 



3 A Green function also has a factor h 1 ^ 2 for each of its external legs, due to the conversions tjj — > i/> and 
A — > A of the corresponding operators. Thus the free electron and photon propagators are of O [Ji] . 
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Hence the propagators can contribute negative powers of a. For diagram (a) in Fig. 4 we 
have A(a) ~ a/q 2 oc I /a. In diagram (b) the four vertices and the dk G d?k of the loop 
integral measure give a 2 - a 2 - a 3 = a 7 in the numerator. The fermion propagators are each 
off-shell by the binding energy and hence oc a~ 2 , while each photon propagator oc a~ 2 as 
in A(d). Thus the four propagators in (b) contribute oc a -8 , implying that A(b) oc a" 1 is 
of the same order in a and h as A{a). This also explains why the ladder sum can diverge at 
the bound state energies, however small is a. The ladder diagrams are enhanced (oc 1/a) 
in the weak coupling limit, and thus give the leading contribution to the bound state poles. 
In diagrams with loop corrections to vertices or propagators the loop momentum does not 
decrease with a, hence such diagrams give true higher order corrections in a and K. 

As a —7- the sum of ladder diagrams in Fig. 4 describes non-relativistic scattering from 
a static Coulomb potential as shown in Fig. 3. Hence the Schrodinger equation gives the 
Born contribution, of lowest order in h, to atomic bound states. The Dirac equation with 
a Coulomb potential A = —a/r can also be obtained from the sum of ladder diagrams 
(with crossed ladders included) in the limit where one of the fermion masses in Fig. 4 
grows large [24] . The relativistic features which distinguish the Dirac from the Schrodinger 
bound states are of the same O (a) as loop contributions. Thus the Dirac equation for 
the Hydrogen atom contains only a subset of the relativistic corrections to the Schrodinger 
equation [20]. 

8. A from Gauss' law 

Having seen that the Schrodinger equation defines a Born term for non-relativistic atoms 
I now return to the hamiltonian field theory formulation (5.1) of bound states. Previously 
I assumed a fixed external potential, now we shall recall how the Coulomb potential arises 
from the equations of motion. I take the constituents to have distinct flavor, an electron 
e~ and a muon [i + for definiteness 4 . 

In Coulomb gauge V • A = the QED lagrangian 

C = —F^F^ + ~ ^ ~ (8-1) 

defines the operator constraint 

-V 2 A°(t, x) = eY J 1>\(t, xWffr x ) ( 8 - 2 ) 

Since there is no time derivative OqA we may express A in terms of the fermion operators 
at each instant of time, 

A°(t,x) = J rf 3 y e _ ^ f (t,y)M^y) (8-3) 

X f=e,H 

4 The annihilation contributions of an e + e~ bound state are in any case of higher order in a than I shall 
be concerned with here. 
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This allows A to be eliminated from the action. From —F^F^/A we get 

- J d 3 x(VA ) 2 = l -J d 3 xA°(-V 2 A°) = ljd 3 x Y^^eA^f (8.4) 

which cancels half the Coulomb interaction of the fermions. The lagrangian becomes 

C = l -{d Q A) 2 - ^FijF» + ~ h^ A ° + ej ■ A - m/)^/ (8.5) 

f=e,n 

with ,4° given by (8.3). So far we made no approximations. 

The e~/i + bound state may be expressed in analogy to the single particle state (5.5) 

as 

\E, t)= I d 3 yi d 3 y 2 4(t, Vx)x{Vx ~ V2)M t > 2/ 2 )|0> (8.6) 

where the rest frame wave function x depends only on the difference of the fermion coordi- 
nates. For small a = e 2 /4-7r the state is non-relativistic and the 4x4 matrix \ is dominated 
by its large ("upper") components. Consequently the couplings of the vector components 
A of the photon field to the fermions are suppressed, and we may at lowest order in a 
neglect them in the hamiltonian, 

H(t) = J d 3 x ^/(*» x ) [ ~ iV • 7 + m + \e^A°{x)] ifj f (t, x) (8.7) 

/=e,/U 

The non-local four-fermion interaction term which arises with the expression (8.3) for A 
looks forbidding, but is harmless in the non-relativistic case where particle production is 
neglected. The only contribution to the eigenvalue equation (5.1) then arises from ip e in the 
hamiltonian annihilating ijr e in the state (8.6) through {i^ a (t, x),tpL(t, y)} = 5 3 (x — y) 5 a p, 
and similarly for the muon fields. Keeping only the two-particle contribution the eigenvalue 
condition becomes 

H(t)\E,t) = J d 3 yi d 3 y 2 ^ e {t, yi ) 7°(-»^l-7 + me)x(Vi-V2) 



x(Vi - V2h°( i V2 + m tA 



a 



m - V2\ 



Requiring the coefficient of ipl(t, yiji^^t, y 2 ) to agree with that of E\E,t) implies 



where 



• 7 + m e ) X (y) - X(yh°(-i^ ■ 7 + ma) = [E - V(\y\)] X (y) (8.9) 



my|) = -j^ (8.io) 

is the standard Coulomb potential. In the weak coupling limit the bound state equation 
(8.9) reduces [8] to the Schrodinger equation 
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In this section I recalled that the ^4° potential can be expressed in terms of the charged 
fields as in (8.3). I next discuss how this expression may be modified by a homogeneous 
solution of the field equations, corresponding to a different boundary condition. This gives 
a linear potential and thus relativistic bound states even at lowest order in a. 

9. A confining boundary condition on A 

The strong interactions are characterized by a scale Aqcd ~ 200 MeV which is not present 
in the classical QCD lagrangian. The scale is introduced when loop corrections are renor- 
malized, resulting in a "dimensional transmutation" where the fixed dimensionless coupling 
is replaced by the running a s (Q). Aqcd also determines the confinement scale through 
vacuum expectation values such as (Q\a s F tlv F tI ' l '\Q) oc Aq CD . It is worth recalling that even 
the QED vacuum is very complicated since the hamiltonian creates and destroys particles, 
implying an infinite number of balancing fluctuations. Nevertheless, in QED one success- 
fully expands around the empty, perturbative vacuum and includes vacuum fluctuations 
only perturbatively. 

The QCD scale is introduced phenomenologically in the quark model, by postulating 
a linear potential with scale k in (2.3). Can there be an approach within perturbative 
QCD that allows to arrive at a quark model picture of hadrons? As I discussed in Sec. 7 
the Schrodinger equation represents a Born level approximation, i.e., it does not include 
loop effects. How could the scale Aqcd appear at Born level in QCD? One possibility 
(and apparently the only one) is via a boundary condition on the color field strength. 
Gauss' law (8.2) (or rather its QCD equivalent) is particularly relevant since it concerns 
the instantaneous field ^4°. Interactions transmitted by the propagating gluon fields A 
necessarily imply that hadrons have gluons in their Fock states, rather than being (at Born 
level) qq and qqq states as in the quark model. The fact that the spectrum of heavy 
quarkonia is qualitatively similar to the atomic spectrum furthermore suggests that they 
bind similarly, through the Coulomb potential A . 

Consider adding a homogeneous solution of Gauss' constraint (8.2) to the standard 
expression (8.3) (still using QED as illustration), 

A°(t,x) = j d 3 y 4 , * , £ ^\{t,y)i, f {t,y) + A 2 l-x (9.1) 

where A and the unit vector £ are independent of x. This contribution implies a con- 
stant energy density ^(VA ) 2 = ^A 2 at asymptotic x. The possibility to add this term 
means that it preserves the stationarity of the action under local variations of ^4°, but 
not necessarily stationarity under the global variation of I. Keeping only the two-particle 
|yi,l/2) — ^(t, 2/i)Vv(*> 2/2)!^) Fock state contribution, 



A°(t,x)\ yi ,y 2 ) 



e ( 1 1 

+ A z £ ■ x 



4vr \\x - i/i| \x - y 2 1 



12/1,2/2) (9-2) 



The electron and muon generate an electric dipole field, which results in the standard 1 jr 
potential in (8.9) (when the infinite, 2/1, 2/2 -independent contributions are discarded). The 
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instantaneous A°-field is thus specific for each position of the fermions. The stationarity of 
the action under the global variation of £ should similarly be ensured separately for each 
Fock state. The field strength contribution may be evaluated as in (8.4) through partial 
integration (except for the A 4 term). This gives to O (e), 



12/1.2/2) 



l - j d 3 x [A A -e^iJ;\(x)A 2 £-x^f(x)\\y 1 ,y 2 ) = ^[A A J d 3 x - eA 2 £ ■ ( Vl - y 2 ) 

f (9-3) 
The first term on the rhs. is proportional to the volume of space and is due to the constant 
energy density ^(VA ) 2 . In bound state calculations this infinite term can be discarded 
provided A is independent of y 1 and y 2 . Stationarity of the second term wrt. variations of 
£ imposes (up to a sign) 

i = (9.4) 
IV1-V2I 

which gives a linear potential contribution to the interaction energy in (8.7), 

V(y 1 ,y 2 ) = 1 2 eA 2 \y 1 -y 2 \ (9.5) 

Assuming a non-vanishing energy density (A 7^ 0) thus leads to a linear potential. It 
is interesting to note that this is consistent with translation invariance only for neutral 
states. For an electron with charge e\ and a muon with charge e% we would have obtained 
a potential V oc \e-iVi ~ e iV2\i which is not invariant under a translation y i — > y i c. The 
restriction to neutral bound states is consistent with color confinement in QCD. 

The linear potential (9.5) is of O (e) and thus leading wrt. the O (e 2 ) photon exchange 
contribution (8.10). This allows to study bound states self-consistently using only the 
linear potential, and introducing photon exchange perturbatively. For relativistic states 
also the propagating components A of the photon field would contribute at O (e 2 ). As I 
mentioned in Sec. 6 the longitudinal components of the photon field must be taken into 
account even for the Hydrogen atom when it is in relativistic CM motion. 

In the previous discussion I considered only two-particle Fock states, as appropriate for 
non- relativistic e~fi + bound states. This restriction can be removed using the observation 
in Sec. 5 that pair production for Dirac states is "hidden" when using retarded boundary 
conditions at asymptotic times. At the Born level, i.e., at lowest order in h, there are no 
loops and the energy eigenvalues are unaffected by the ie prescription of the propagators. It 
should be recalled that despite the apparent simplicity of the ensuing, Dirac-type "valence" 
wave functions the relativistic bound states actually contain an infinite number of pairs 
when Feynman boundary conditions are imposed at t = ±00. 

The generalization of the "retarded vacuum" (5.3) to the case of two fermions e~ and 
li + is straightforward, 

\0) R = N- 1 H4(p,\)bl(p,\)\0) (9.6) 
P ,x 

implying 

Mx)\Q)r = ^1(x)\0) r = (9.7) 

When the \E,t) bound state in (8.6) is defined using the retarded \0)r rather than the 
perturbative vacuum |0) the property (9.7) ensures that only two-particle states contribute. 
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10. qq states in QCD 

I used QED in the above discussion to illustrate the essential features of the approach. The 
method of deriving bound states at Born level in h and at lowest order in the coupling can 
also be applied to QCD, with 

Cqcd = -\F^F* U + ^j{i$-gA a T^ B -m } )^f 

f=u,d 

Fr = &>Al - d v A» - gf abc A%A» (10.1) 
The retarded vacuum now includes a product over color to make it gauge invariant, 

\0)r = N- 1 J] ^t(p,A)6f(p,A)|0) (10.2) 

p,X,A 

and the ud meson state is 

\E,t) =J2 f d 3 y 1 d 3 y 2 ^(t,y 1 )x AB (y 1 ,y 2 )^(t,y2)\0)R (10-3) 

A,B ^ 

This state is invariant under (time-independent) gauge transformations ip(t, x) — > U(x)ip(t, x) 
provided that the wave function is transformed as 

x(vi,y 2 ) -> U(y 1 ) x (y 1 ,y 2 )U\y 2 ) (10.4) 

We may then look for stationary states of the form (10.3) in a Coulomb gauge where the 
wave function has the standard "color singlet" form, 

x AB (yi,y2) = s AB x(yi,y2) (io.5) 

The equation of motion for the gluon field 

^QCD = {d»8 ab + gfa b cA^F b w - g £ ffipT^f = (10.6) 

has the perturbative solution 

Al(t,x)=A%-x+ f d 3 y g _ ^f(t, y)T^ B ^f(t,y) + O (g 2 ) 
J n\ V\ f 

A a (t,x) = 0(g) (10.7) 

provided A a = for a 7^ 3,8. Introducing a homogeneous solution only for the color 
diagonal generators T3 and T$ is necessary to ensure the color structure (10.5) of the wave 
function and implies that the non-abelian terms in (10.6) are of O (g 2 )- 

For the solution (10.7) the interaction part of the Lagrange function Lj(t) = f d 3 x Cj^qcd 

is 

Lj(t) = ljd 3 xJ2 [K - gA 2 a h ■ x ]T tf\t, x)T^ B i;f(t, x) + O {g 2 ) ] (10.8) 
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The (infinite) energy proportional to the volume of space must be the same for all Fock 
components of the bound state (10.3), which requires that the field strength 

A 4 - £ A * ( 10 - 9 ) 

a=3,8 

is a universal constant. I do not show this term in the following. The remaining part of Lj(t) 
in (10.8) should be stationary on each Fock component I2/1, 2/2; C) = V'^(2/i)V ; rf'(2/2)|0)-R> 

Li(t)\ yi ,y 2 ; C) = -f £ A* if c 4 ■ (y x - y 2 )\y lt y 2 ; C)+0 [g 2 ) (10.10) 

a 

The variation wrt. the direction of £ a imposes (up to a sign), 

I _ T a C {V\-V2) no 111 

The Lagrange function (10.10) must be stationary also wrt. variations of the ratio A|/A|, 
under the constraint that A in (10.9) is fixed. Together with (10.11) this gives 

Lj(t)\ yi ,y 2] C) = - 9 ^J2 a (T^) 2 \ Vl - y 2 \\ yi ,y 2 ;C) + 0(g 2 ) (10.12) 

The SU(3) generator identity 

E a T a AB T^ D = \{5 AD 5 BC - \8ab8cd) (10.13) 

gives ^ a (Ta C ) 2 = \i ensuring that the eigenvalue in (10.12) is independent of the color 
C of the quarks. For the interaction hamiltonian Hj(t) = —Lj{t) + O (g 2 ) we have 

0A 2 

Hj(t)\ yi ,y 2 ;C) = ^j=\vi -V 2 \ \Vi,V*C) + O {g 2 ) (10.14) 

The fact that the linear potential has the same strength for all color components of the 
wave function is consistent with the diagonal ansatz (10.5). 

The eigenvalue condition for the bound states defined by (10.3) and (10.5) is then, up 
to terms of O (g 2 ) , 

H(t)\E,t) = [ d 3 xd 3 yi d 3 y 2 V ^f(t,x)\-iV x -j + m f + ^=\( yi - y2 )\ 7 °] itf(t,x) 

J f,A,B 1 2 ^ 3 J ' 

x ^(t,y 1 ) x (y 1 ,y 2 )^(t,y 2 )\0)R = E\E,t) (10.15) 

Using 

{ifif^t, x),^(t, y)} = S 3 (x - y) 5 aP 6 ff , 5 AB and if£(x)\0) R = 4>f(x)\0) R = 

(10.16) 

we get the condition for the wave function x{yx^y 2 ) OI " (10-5), 

7 (-^ r 7+m u )x(yi,y 2 )-x(y 1 ,2/ 2 )7 (^2+m d ) = [E-V(\ yi -y 2 \)}x(yi,y 2 ) (10.17) 
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with the linear potential 5 

V(\y\) = ^=\y\ + 0(g 2 ) (10.18) 

This bound state condition for the wave function is valid in any frame. For a state 
(10.3) with CM momentum k the wave function has the form 

*(Vi, V 2 ) = e^ + y^ 2 Xk ( yi - y 2 ) (10.19) 

Substituting this in (10.17) gives 

-iV-[cx,Xk(y)} + lk-{<x,Xk{y)}+m ul x k (y)-Xk{yh°m d = [E k -V{\y\)] Xk (y) (10.20) 

where a = 7°7- As I discussed in Sec. 6, equal-time wave functions transform dynamically 
under boosts. Thus (10.20) is not manifestly Lorentz covariant. Nevertheless, since we 
arrived at this equation through an expansion in the fundamental parameters h and a s , 
the energy eigenvalues should have the correct dependence on the CM momentum: E k = 
\/ k 2 + M 2 . Remarkably, this relation is found to hold [25]. The wave function transforms 
dynamically under boosts, Lorentz contracting at a rate oc 1/(E — V) which depends on 
the potential, and hence on the separation \y\ between the quarks. 

The bound state equation (10.20) is (to my knowledge) the only case where a correct 
fc-dependence of the energy has been obtained. This occurs only in the case of an exactly 
linear potential, as expected since the O (g 2 ) gluon exchange corrections were neglected in 
the derivation. 

11. qqq states in QCD 

The derivation of uds bound states (baryons) is similar to the above qq case (mesons), but 
contains some new elements. The state is defined as 

\E,t)= ^ f d'y^y^y^^y^^ 
A,B.C J 

(11.1) 

where now 

|0>* = N- 1 H d^(p, A) df (p, A) df (p, A)|0) (11.2) 
p,A,A 

The baryon state (11-1) is invariant under time independent gauge transformations tp A (t, x) 
U AA ' (x)tfj A ' (t, x) provided the wave function is transformed as 

Xabc(xi,x 2 ,x 3 ) -> U AA ' ( Xl )U BB ' (x 2 )U cc ' (x 3 )xa'B'C'(xi,x 2 ,x 3 ) (11.3) 

I assume that there is a gauge where 

x7Sh a3 {vi^yz) = e A Bcx aia2a:i (yi,y2,y 3 ) (n-4) 

The coefficient differs from that of Eq. (5.25) in [8] due to a different regularization of the infinite 
contribution oc the volume of space in (10.8). 
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The solution (10.7) of the QCD equations of motion gives the Lagrange function (10.8) 
with the same divergent contribution as before which requires that A in (10.9) is universal 
(independent of the state). Applying the Lagrange function (minus the divergent term) to 
a specific Fock component of the bound state (11.1) where, according to (11.4) A,B,C is 
some permutation of 1,2,3, gives 

L(t)\ yi ,y 2 ,y 3 ;ABC) = L(t)^ (y,)^ (y 2 )^ (y 3 )\0) R (11.5) 

= - § E A * *° ■ ( T « A yi + T * B yi + T - C yz) \v^v* v* ABC ) 

a 

The condition of stationarity aligns t a with the vector it multiplies in (11.5), hence 
L(t)\ yi ,y 2 ,y 3 ;ABC) = -^^T^yx+T^y^ + T^y^y^-ABC) (11.6) 

a 

Under translations y i — > + c the vector Y a = T l AA y 1 + T BB y 2 + T^ c y 3 shifts by 
(Tr T a ) c = 0. Hence the color structure (11.4) ensures translation invariance. 

The expression (11.6) is extremal for A|/A| = | "5:^3 1 / 1 Ir^g | , which depends on the quark 
positions y i . The interaction hamiltonian, i.e., the potential, is then 

2 



V(y 1 ,y 2 ,y 3 ) = 9 -^ r ^Y 2 3 + Yl (11.7) 

In order to be consistent with the color structure (11.4) this potential should be independent 
of the specific color permutation A, B, C of the quarks at position y\,y 2 ,y 3 - This may be 
seen as follows: 

£ Y l = £ [ T a A yi + T a B V2 ~ (T a AA + T? B )y 3 ] 2 (11.8) 

a a 

= E [( T a A ) 2 (yi - vsf + (T a BB ) 2 (v2 - y,) 2 + 2T a AA T a BB ( yi - y a ) • (y a - Vs\ 

a 

where according to (10.13) { T a A f = \ and Y.a 2T a AT a B = ~l (for A ^ B). Hence 
the baryon potential is independent of the color permutation A,B,C and equals 

V( yi ,y 2 ,y 3 ) = 0=V(Vi ~ + (» 2 - Vs) 2 + Gfo ~ l/i) 2 (H-9) 

In the limit where two quarks are in the same position, e.g., y 2 = y 3 , this potential 
coincides with the meson potential (10.18). An analogous potential cannot be constructed 
for more than three quarks, as there are then more independent quark separations than 
diagonal SU(3) generators. 

Proceeding as in the qq case results in the bound state condition on the wave function 
(11.4) [8], 

3 

E [-A-iVj "1, + mj)] x(Wi,l/2,V 3 ) = (E - V) X {Vi,V 2 ,Vz) (U-10) 

i=i 

with V given by (11.9). 
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12. Discussion 



The above presentation covered much ground. Some of the topics were of a very basic and 
established nature, while other parts may be classified as recent research proposals. To 
make it easier to digest the arguments I recapitulate and comment on the main steps. 

I discussed a perturbative approach to bound states in field theory. For QED atoms we 
know that this is the right method, since it gives accurate agreement with data. Applying 
perturbation theory to bound states requires somewhat different methods than used for 
scattering amplitude calculations. In order to generate a bound state pole one needs to 
sum an infinite set of Feynman diagrams. In the a — > limit ladder diagrams of the type 
shown in Fig. 4 give the leading contribution when the external momenta are scaled with 
a so as to be compatible with the bound state wave function. The residue of the pole 
satisfies the Schrodinger equation, which defines the Born term in a systematic evaluation 
of higher order corrections [19, 26]. 

Perturbation theory can give us insights into the properties of relativistic bound states 
from "first principles" , and its relevance for QCD bound states merits careful consideration. 
There are indications that the QCD coupling has a moderate value a® ~ 0.5 even at 
distances of O (1 fm). Many issues involving relativistic states have received less attention 
than they deserve, such as the boost dependence of the Hydrogen atom wave function in 
QED [7]. As I noted in Sec. 6, this offers an example of how the Lorentz contraction of a 
classical stick is realized, and differs, in quantum mechanics. 

In perturbation theory already the lowest order contribution is expected to provide a 
good qualitative description of the exact result. Hence I focussed on the physical principles 
of bound states at Born level. The Hydrogen atom may be thought of as bound by the 
instantaneous Coulomb field A (x; x i, X2 ) which is determined by Gauss' law (8.2) sepa- 
rately and differently for each position Xi, x<i of the electron and proton. Both the electron 
and the proton interact with this Coulomb field - double counting is prevented because the 
field energy (8.4) cancels half of the interaction energy. The Coulomb field A measured 
by an external observer far from the bound state is then given by the coherent sum over all 
configurations of the atom, A°(x) = J dx\dx2 i p[x\^ X2)A°(x; x±, x<i) where tp is the wave 
function. 

In courses on relativistic quantum mechanics the Dirac equation is often introduced 
through its historical context, with a wave function <p(x) that is a c-number. In a field 
theory context the Dirac equation later reappears as an exact operator-valued equation 
derived from the QED action. The similar notation obscures the fact that the c-numbered 
Dirac equation is only an approximate relation which can be derived from summing a 
subset of Feynman diagrams in a certain limit [24]. The negative energy components of 
the Dirac wave function <p(x) are often vaguely explained as related to antiparticle effects. 
In Sees. 4 and 5 we saw that <p(x) is the wave function obtained using retarded propagators 
and describes a state which, with Feynman propagators, has contributions from an infinite 
number of particle-antiparticle pairs. The two descriptions give exactly the same bound 
state energy in the absence of loop corrections. This possibility to describe a state with 
infinitely many pairs using a single-particle wave function might shed light on the puzzle 
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of why hadron quantum numbers reflect their valence quark degrees of freedom only, with 
no reference to the sea quarks and gluons. 

The use of retarded propagators is only possible in the absence of loops, which in turn is 
generally understood to imply Born level, or lowest order in K. The concept of a Born term 
for bound states appears not to have been discussed before [8, 9]. The Schrodinger equation 
is the Born term for non-relativistic states (Sec. 7). The h expansion is equally applicable 
to relativistic bound states, for which it provides a well-defined first approximation. Since 
H is a fundamental parameter each order in the expansion must have all the symmetries of 
the exact result. As I discussed in Sec. 6, the full Poincare invariance of quantum theory is 
not manifest due to the choice of a quantization surface {e.g., equal time or equal light-front 
time) [21]. Hence the Born term of equal-time bound states has a hidden, dynamical boost 
invariance. This can hardly be obtained without a precise theoretical framework such as 
perturbation theory [27]. 

Students should be aware that my proposal in Sec. 9 to consider a homogeneous so- 
lution to Gauss' law (Eq. (9.1)) is an "educated speculation". There is not sufficient 
experience yet to tell whether such a solution gives fully self-consistent and meaningful 
results. It fulfills the essential requirement of giving a stationary action, thus maintaining 
the symmetries of the theory. The bound state equation has manifest rotational symmetry 
in the CM and dynamic boost invariance for the purely linear potential of the Born term. 
The application to qq and qqq bound states in QCD is straightforward and maintains 
color covariance. Further work is needed to study the properties of such solutions. The 
present hamiltonian formulation should allow to study many physical observables such as 
form factors, parton distributions and scattering phenomena. The relevance of imposing 
the boundary condition (10.7) on A°{x) which gives a finite energy density as \x\ — > oo 
will depend on the outcome of those applications. 

I should emphasize that the present approach does not explain why there is confinement 
in QCD but not in QED. The vacuum energy density A 4 in (10.9) is a free parameter, which 
must be set to zero "by hand" for QED. This is only a proposal for how confinement can 
be described in a field theory context. 

I did not discuss the important issue of chiral symmetry. In the spirit of the present 
approach a possibility to describe spontaneous symmetry breaking may be to impose chi- 
rally non-symmetric states as a boundary condition, in analogy to the non- vanishing of the 
energy density, A ^ in (10.9). 
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